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PREFACE 


This  report  was  prepared  for  the  U.S.  Army  Night  Vision  ar.d 
Electro-Optics  Laboratory,  Laser  Division,  Fort  Bolvoir,,  Virginia, 
it  constitutes  1’i'ci fic-Sierra  Research  Corporation's  final  report  on 
contract  DAAK20-79-C-0040  with  that  agency.  An  earlier  version  was 
presented  as  a  paper  to  the  Fourth  Smoke /Obscurants  Symposium,  held 
22-23  April  .1980  at  Harry  Diamond  Laboratories,  Adolphi,  Maryland. 


SUMMARY 


This  report  develops  an  algorithm  for  inferring  the  absolute 
transmittance  and  extinction  coefficient  throughout  an  obscurant 
cloud  from  which  a  resolvable  signal  can  be  detected,  using  only  rela¬ 
tive  ambient  and  cloud  lidar  signatures.  Simple  expressions  are 
derived  to  approximate  the  errors  in  the  inferred  values  for  given 
uncertainties  of  the  cloud  parameter.  The  obscurant  mass  concentra¬ 
tion  and  concentration  path  integral  (CI,  values)  are  shown  to  be 
proportional  to  the  extinction  coefficient,  so  that  relative  values 
of  mass  concentration  and  CL  values  can  be  directly  inferred  from 
the  extinction  coefficient.  A  point  calibration  measurement  is  needed 
to  infer  absolute  values,  however. 

The  results  are  applied  to  actual  field  data  for  dust  clouds 
measured  during  the  DIRT-I  experiment. 
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I.  INTRODUCTION 


Work  perforat'd  under  a  previous  contract  [Warren  and  Lutomirski, 
1979]  established  the  validity  of  an  algorithm  that  inverts  lidar 
signatures  for  use  in  characterizing  the  aerosol  dynamics  of  obscurant 
clouds  under  battlefield  conditions.  The  algorithm  integrates  a  time- 
resolved  lidar  signature  to  estimate  (.1)  transmittance,  (2)  extinction 
coefficient,  (3)  mass  concentration,  and  (4)  concentration  path  inte¬ 
gral  (CL  value)  for  spatially  varying  obscurant  clouds.  In  the  earlier 
study,  computer  simulations  of  the  lidar  return  used  to  test  the  al¬ 
gorithm's  sensitivity  to  noise  and  parameter  uncertainties  (such  as  in 
the  backseat ter-tc-extinction  ratio)  showed  that  locally  measuring  the 
scattering  properties  of  the  obscurant  cloud  would  allow  the  mass  con¬ 
centration  and  CL  values  to  he  inferred  throughout  the.  obscurant  cloud 
where  the  lidar  return  is  above  the  noise. 

The  present  effort  establishes  the  feasibility  of  inferring  the 
transmittance  and  the  extinction  coefficient  throughout  the  cloud  using 
only  the  ambient  and  obscurant  cloud  lidar  signatures,  without  addi¬ 
tional  point  calibration  measurements.  We  find  that  the  required  nor¬ 
malization  can  be  accomplished  by  vising  the  ambient  signature  and 
either  waiting  until  the  cloud  is  tenuous  or  having  a  scanning  capa¬ 
bility  to  probe  the  region  behind  the  cloud.  A  target  behind  the  cloud 
is  useful  for  enhancing  the  reliability  of  the  normalization,  but  is 
not  essential.  The  relative  mass  concentration  and  CL  value  of  an 
obscurant  cloud  can  be  inferred  directly  from  the  extinction  coef¬ 
ficient.  The  absolute  mass  concentration  and  CL  value  can  be  found, 
if  needed,  by  performing  <?.  single  point  measurement  along  the  lidar 
path  within  the  cloud.  Directly  measuring  the  concentration  instead 
of  measuring  the  particle  size  distribution  is  preferable  for  that 
purpose . 

Applying  the  above  method  to  lidar  signatures  from  dust  clouds 
measured  during  the  DIRT- I  tests  [van  der  Laan,  1979],  sample  calcu¬ 
lations  show  maximum  errors  of  less  than  5  percent  in  transmittance 
and  less  than  20  percent  in  the  extinction  coefficient  for  a  10  percent 


v  ~  )  V 


uncertainty  in  the  inferred  scattering  parameters.  For  the  same  un¬ 
certainty  in  the  scattering  parameters,  the  errors  in  the  concentra¬ 
tion  and  the  CL  value  are  less  than  20  percent  and  12  percent,  re¬ 
spectively.  The  report  establishes  that  simple  equations  can  predict 
these  errors. 

Section  II  defines  the  assumptions  for  the  model  and  derives 
expressions  for  the  transmittance  and  extinction  of  the  cloud,  giving 
solution  methods  with  or  without  a  point  measurement.  Section  III 
investigates  the  sensitivity  of  the  algorithm  to  uncertainties  in  the 
obscurant  cloud  backscatter-to-extinction  ratio  and  derives  simple 
closed-form  estimates  for  the  uncertainty  in  the  estimated  transmit¬ 
tance  and  extinction  coefficient. 

Section  IV  calculates  the  concentration  and  CL  values  of  the 
smoke  or  dust  (hereafter  simply  smoke)  and  recommends  an  alternative 
to  obtaining  local  measurements  of  the  particle  number  distribution. 
Section  V  applies  the  inversion  algorithm  to  data  taken  during  DIRT-1, 
numerically  calculating  the  transmittance,  extinction  coefficient,  mass 
concentration,  and  CL  values  and  estimating  the  ranges  of  possible 
errors.  Section  VI  summarizes  the  proposed  model  and  the  results  of 
applying  the  algorithm  to  the  DIRT-I  data. 

The  need  for  a  practical  approach  to  deriving  detailed  informa¬ 
tion  on  the  spatial  and  temporal  dynamics  of  battlefield  obscurants 
is  well  established.  An  additional  strong  motivation  for  the  analysis 
described  here  is  the  concurrent  development  by  the  Army  Night  Vision 
and  Electro-Optics  Laboratory  of  a  two-color,  three-dimensional  scan¬ 
ning  lidar.  The  added  flexibility  of  that  system  substantially  en¬ 
hances  the  usefulness  of  the  inversion  algorithm  in  that  the  required 
trar smittance  normalization  can  be  performed  by  selecting  an  optimal 
combination  of  spatial  scan,  wavelength,  and  temporal  signature.  The 
work  here  represents  a  step  toward  exploiting  this  improved  capability 
for  inferring  smoke/dust  dynamics  in  the  field. 
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II.  ALGORITHM  FOR  ESTIMATING  OBSCURANT 
TRANSMITTANCE  AND  EXTINCTION 


The  formal  solution  to  the  radiative  transport  equation  for  a 
spatially  varying  medium  was  solved  for  a  general  scattering  medium 
by  Warren  and  I.utomlrski  [1979].  Their  solution  can  be  specialized 
to  lidars,  resulting  in  the  single-scattering  lidar  equation 


where.  P(z)  is  the  instantaneous  backscatter  power  received  at  time  t 
[proportional  to  the  range  z(t  =  2z/c)]  and  A  is  a  constant  for  given 
lidar  parameters — which  include  wavelength  of  the.  lidar  laser,  losses 
or  efficiency  of  the  transmitting  and  receiving  optics,  effective  re¬ 
ceiver  aperture,  backscatter  phase  function,  and  a  dimensional  con¬ 
stant.  The  backscatter  coefficient  at  range  z  is  given  by  3(z);  the 

one-way  transmittance  of  the  medium  to  range  z  is  given  by  T(z).  The 

2 

two-way,  round-trip  transmittance  of  the  medium  to  range  z,  T  (z) ,  can 
be  written  as 


where  c(z)  is  the  extinction  coefficient  at  range  z.  The  center-to- 
center  separation  of  the  transmitter  and  receiver  is  b,  and 

C(z)  =  a2  +  02z2  +  (aT  +  0 ,f z ) 2  ,  (3) 


where  a^  is  the  receiver  aperture  radius,  0R  is  the  f ield-of-view  half- 
angle,  a,^  is  the  transmitter  aperture  radius,  and  0^  is  i  lie  trans¬ 
mitter  beam  divergence  half-angle.  Equation  (3)  assumes  a  Gaussian 
aperture  and  a  Gaussian  f ield-of-view  receiver. 


> t,uu 
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Equation  (1),  the  lidar  equation,  is  valid  as  long  as  the  trans¬ 
mitter  has  a  short  pulse  (less  than  a  few  hundred  nanoseconds)  and  a 
narrow  beam  width,  (at  most  a  few  centimeters);  a  small  divergence  (on 
the  order  of  a  few  milliradians)  is  required  for  there  to  be  no  sig¬ 
nificant  transverse  variations  in  the  lidar  backscatter  power  return. 
The  receiver  f ield-of-vicw  must  be  narrow  (on  the  order  of  a  few 
milliradians)  to  minimize  the  effects  of  multiple  scattering. 

The  lidar  backscatter  power  return  from  a  smoke  cloud  raised  by 
an  explosion  is  composed  of  two  parts — the  ambient  air  and  the  smoke 
cloud  return.  The  cloud  usually  comprises  a  range  of  particle  sizes 
describable  by  a  particle  size  distribution  n(r,  z) ,  where  r  is  the 
radius  of  the  particle.  In  what  follows,  the  subscript  a  refers  to 
ambient  air,  s  to  smoke  particles  only,  and  c  to  ambient  air  and  smoke 
combined . 

To  use  Eq.  (1),  we  make  several  assumptions  about  the  backscatter 

and  extinction  coefficients.  First,  we  assume  that  the  ambient  air 

backscatter  coefficient  3,  and  extinction  coefficient  e  are  constant 

d  d 

throughout  the  range  and  are  unaffected  by  the  smoke.  In  other  words, 
the  partial  volume  occupied  by  the  smoke  is  small  and  the  relaxation 
time  for  3.  due  to  the  explosion  of  the  smoke  is  fairly  short,  so  that 

d 

3a  and  return  to  their  equilibrium  values  soon  after  the  explosion. 

Second,  we  assume  that  for  the  ambient  air  where  the  coefficients 
are  constant,  3  “at.  As  shown  in  Fenn  [1966],  the  relationship 

Q  3d 

between  the  extinction  and  backscatter  coefficients  is 


3  “  ae  ,  (A) 

where  both  a  and  l  are  constants.  Since  the  lidar  equation  [Eq.  (1)] 

is  for  tlie  single-scattering  model,  l  is  taken  to  be.  1,  giving  a  linear 

relationship  between  the  backscatter  and  the  extinction  coefficients. 

Finally,  we  assume  that  n  (r,  z,  t) — the  number  of  smoke  particles 

s 

with  radius  r,  range  z,  and  time  t — is  separable,  such  that 


n  (r,  z,  t)  “  n.  (r)n,  (z,  t)  , 

;•»  .1.  S  4.  %  ? 
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and  that  the  particle  size  distribution  is  constant.  The  total  number 

of  particles  per  unit  volume  n,,s(z,  t)  may  vary  with  range  and  time, 

which  will  cause  £  (z)  to  so  vary.  That  assumption  allows  expressing 
s 

the  smoke  backseatter  coefficient  as  B  (z)  =  a  e  (z) ,  where  a  is 

s  s  s  s 

independent  of  range  and  time. 

The  lidar  backseatter  power  return  in  ambient  air  is 


A3  T?(z) 

■n  /  \  3.  cl 

Vz>  ’  ow"  “tr> 


(* 1 v 

\G(z)  )  ’ 


(6) 


in  the  presence  of  smoke,  the  total  backseatter  power  return  becomes 


AB  (z)T"(z) 

V2)  =  — tloo - exp 


m- 


(7) 


The  backseatter  and  extinction  coefficients  for  combined  ambient  air 
and  smoke  consist  of  the  sum  of  the  component  coefficients: 


4 

•  * ! 
i  - 


i 


Bc(z)  =  Ba  +  0s(z)  ;  ec(z)  “  ea  +  Eg(z)  .  (8) 


Combining  Eqs.  (7)  and  (8)  results  in  the  following  relation  for 

T^z): 


T^(z)  =  Tj(z)T^(z)  .  (9) 


Dishing  Kq.  (7)  by  Eq.  (6)  and  substituting  the  expressions  for 
Bc  and  T^  gives 


p  U)  2  3  (Z) 

pfuT  ■  V*>  +  S-  i (z)  • 

3  3 


(10) 
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Now  that  no  explicit  knowledge  of  the  actual  ltdnr  parameters  (such 
as  wavelength ,  transmitter  or  receiver  aperture,  efficiency  of  the 
lidor  pulse,  power)  is  necessary;  all  such  information  is  carried  in 
the  ambient  backscatter  signature  P„.  A  Gaussian  aperture  and  field 
of  view  were  assumed  for  G<k>,  but  since  GOO  cancels,  Eq.  (10)  is  now 
valid  for  more  general  problems-tbot  Is,  those  that  have  the  form  of 

Eq.  (1)*  ,y 

The  transmittance  of  the  smoke  T>)  can  be  differentiated  with 
respect  to  a.  Substituting  6 „<«)*>.  for  S/x)  then  gives  an  expres- 

sion  for  0s(z) : 


B  . - 

2-r(z) 


d'i’2(c) 


substituting  Eg.  (11)  into  Eq.  (10)  gives  a  first-order  differed- 


tie'll  equation  for  Tg: 


V  U)  2  %  /  dTs(z) 

TJzJ  *  Va>  +  2Ba  '  dz 


which  can  be  integrated  to  yield 


P  (a1) 


:Zs(z)  =  [exp  (qz)]  1  “  q  f  ^  ’ 

J  n 


where  q  =  ^a/%  a"d  '^(0)  ''  U  Hange  zer°  “  0)  C5n  bG  ch°Sen 
arbitrarily—Vmt,  for  simplicity,  it  should  be  chosen  so  that  z  =  0 

is  in  a  region  in  front  of  the  smoke  cloud,  i.e.,  where  Pc<z)  -  Vz)  1 
Thus,  the  exact  position  of  the  lidar  is  not  needed  for  the  calcula¬ 
tion  of  T2(s0.  Once  tJ(z)  has  been  calculated,  Eq ,  (10)  can  be  solved 
for  the  backscatter  coefficient;  substituting  for  &s*  the  extinction 
coefficient  at  range  z  is 
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e 

s 


(z) 


a 

2 


(14) 


The  smoke  cloud  data  usually  consist  of  the  range-resolved  ampli¬ 
tude  of  the  backscatter  power  return  taken  at  predetermined  time  inter¬ 
vals,  beginning  before  the  smoke  cloud  is  introduced  and  lasting  until 
it  dissipates.  Figure  1  shows  the  10.6  pm  backscatter  data  from  evert 
F-6  during  DIRT-I.  The  range  (in  meters)  is  given  along  the  abscissa, 
and  the  log  of  the  backscatter  amplitude,  along  the  ordinate.  At 
T  =  0,  the  spike  at  100  m  is  due  to  a  reflector  used  to  calibrate  the 
range;  the  spike  at  2000  m  is  a  target  used  to  calibrate  the  range  and 
measure  transmittance.  The  ambient  air  backscatter  amplitude  P  (z)  is 

cl 

shown  in  the  T  =  0  frame.  In  the  next  frame,  T  +  19s,  detonation  of 
the  explosive  has  already  occurred,  causing  an  extremely  dense  cloud 
at  1000  m;  there  is  therefore  little  penetration  of  the  laser  beam 
into  the  cloud,  and  the  target  at  2000  m  is  completely  obscured.  As 
time,  passes,  there  is  greater  penetration  until  the  cloud  is  tenuous 
enough  to  allow  t he  laser  beam  to  pass  through,  and  the  backscatter 
amplitude  from  the  target  behind  the  cloud  becomes  visible. 

The  smoke  cloud  daU  <e  P  (z)  and  P  (z) ,  but  a  value  for  q  must 
be  obtained  to  solve  for  the  smoke  transmittance  and  the  extinction 
•oefficient.  There  are  two  main  methods.  Method  I  is  to  take  a  local 
point  measurement  of  q  along  the  path  of  propagation  of  the  lidar  beam. 
Method  Ila  is  to  pick  a  set  of  backscatter  amplitude  data  in  which  the 
smoke  is  prevalent  but  the  backscatter  amplitude  of  the  ambient  air 
from  behind  the  cloud  is  both  visible  and  well  above  the  noise  level. 


Far  beyond  the  smoke  cloud,  where  z  =  z  ,  there  is  essentially  no 

smoke;  so  the  extinction  coefficient  for  the  smoke  c  (z  )  is  zero  in 

s  ns 

that  region,  and  the  ratio  P  (z  )/P  (z  )  is  constant.  From  either 
b  c  ns  a  ns  „ 

Ec, .  (L2)  or  Eq.  (14)  it  follows  that  P  (z  )/P  (z  )  equals  T  (z  ) 

1  1  c  ris  a  ns  s  ns 

in  the  region  beyond  the  smoke  cloud.  Titus,  the  value  of  q  can  be 

2  2 

found  by  varying  q  in  Eq.  (13)  until  T  (z)  equals  T,(z  ).  An  alterna- 

g  s  ns 

tive  method  (lib)  is  to  place  a  target  in  the  smoke-free  region  beyond 
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Fig.  1--Event  F-6  10.6  ym  backscatter  data 
from  DIRT-I  [van  dor  Laan,  1979]. 
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thi'  smoke  cloud  in  line  with  the  loser  beam.  The  ratio  of  the  back- 

scatter  amplitude  of  the  target  In  ambient  air  (R  )  to  the  backscatter 

a 

amplitude  of  the  target  with  the  smoke  cloud  wlu  n  th^  target  is  clearly 

2 

visible  (R  )  is  just  T  (z  )  ,  where  z  is  the  target  range,  i.o., 

T„(x.t)  =  Rt,/Ra. 

The  T~  found  either  way  should  result  in  identical  values,  but 
h  2 

method  1Tb  will  Improve  the  accuracy  of  T  due  to  the  higher  signal- 

S  2 

to-noise  ratio  (SNR)  of  the  target  backscatter  data.  The  value  of  T 

s 

can  he.  substituted  into  Kq.  (13)  and  combined  with  P  (z)  and  Pfl(z)  to 

iterate  for  a  value  of  q.  In  Llie  lidar  equation,  P(z)  is  in  absolute 

units  such  as  watts,  but  by  using  the  ratio  of  1’  to  1’^,  only  the 

relative  values  of  P  and  P  are  needl’d  to  .solve  Cor  the  smoke  trans- 

c  a 

mittanee  and  the  extinction  coefficient.  Also  note  that  if  the  value 

for  6  is  known,  methods  I  fa  and  Tib  supply  a  technique  to  find  the 
a 

value  of  remotely. 

If  q  is  constant  for  a  given  smoke  cloud,  then  once  the  value  of 

q  has  been  found,  T  (z)  and  e  (z)  can  be  calculated  for  all  times  over 

s  s 

which  the*  backscatter  data  are  taken.  Tf  the  point  measurement  method 
2 

is  used,  T  ,(z)  and  c  (z)  can  be  calculated  for  very  dense  clouds  u> 

S  H 

to  the  deepest  point  of  penetration  of  the  .laser  beam;  while,  with 

method  ITa  or  Tib,  the  cloud  has  to  dissipate  so  that  backscatter  data 

can  be  collected  from  beyond  it.  If  the  lidar  is  capable  of  scanning 

the  cloud  horizontally  or  vertically,  then  either  method  I  or  method 

I  la  can  find  q  when  the  central  region  of  the  cloud  is  too  dense  for 

tin?  laser  beam  to  penetrate.  By  method  I,  q  can  always  be  found.  A 

scanning  lidar  can  collect  backscatter  data  from  behind  the  less  dense 

edges  of  tin!  cloud  where  there  is  no  smoke.  Then  q  can  be  calculated 

by  method  Tla  and  the  remainder  of  the  cloud  can  be  scanned  to  obtain 

2 

a  three-dimensional  estimate  of  T  (z)  and  £  (z)  for  all  points  from 

s  s 

which  a  signal  above  the  noise  level  Is  received. 


ITT.  SENSITIVITY  ANALYSIS  OK  Till-  TRANSMIT!  ANCK 


AND  THE  EXTINCTION  COEFFICIENT 


The  transmittance  and  the  extinction  coefficient  were  calculated 
in  Eos.  (13)  and  (14).  Both  functions  are  dependent  on  the  range  z , 
hut  errors  In  range  are  usually  minimal.  (Note  that  range  and  time 
resolution  are  related  by  At  =  2Az/c,  so  that  a  10  m  range  resolution 
requires  a  60  ns  time  resolution  in  the  recorded  Ildar  signal.)  Ig¬ 
noring  errors  due  to  range,  the  uncertainty  in  the  transmittance  is 
due  to  noise  in  the  hackscattor  power  returns  of  the  ambient  air  and 

smoke  cloud  (P  and  P  )  and  the  coefficient  a .  In  Warren  and 
a  c 

Lutomirski's  [1979]  work  on  simulated  Ildar  signals  with  superimposed 
Gaussian  noise,  numerical  computer  output  showed  that  the  noise¬ 
averaging  produced  by  the  Integration  for  transmission  gives  reason¬ 
ably  correct  results  for  low  SNRs — for  example,  for  short-range  trans¬ 
mittance  of  signals  having  an  SNR  f:  .1.  The  major  source  of  error  for 
the  transmittance  is  therefore  determining  q.  Differentiating  the 
expression  for  transmittance  with  respect  to  q  gives  an  expression  for 
the  uncertainty: 


AT 
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4-  qz  +  4,'  ( 
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qz 
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(15) 


Equation  (15)  Is  difficult  to  apply,  but  for  regions  where 
qz  «  1, 


At 

S 

— —  —  r  ^ 

T  ‘  ‘ 
8 


(16) 


Figure  2  plots  Kq .  (lb)  for  various  values  of  T^ .  The  estimate  of 
is  not  very  sensitive  to  errors  In  q  for  large  values  of  T^.  For  small 
values  (which  correspond  to  denser  clouds),  T  becomes  much  more  sensi¬ 
tive  to  errors  in  q,  implying  that  for  dense  clouds,  local  measurements 


of  q  may  be  necessary  to  minimize  uncertainty  in  q.  For  example,  in 
a  dense  cloud  where  the  one-way  transmittance  of  the  smoke  is  only 
.25  and  the  acceptable  error  in  is  10  percent,  the  allowable  error 
in  q  is  about  1  percent. 

The  errors  for  the  extinction  coefficient  are  more  serious  than 

those  for  the  transmittance.  As  Eq.  (14)  shows,  there  is  no  integral 

to  average  out  noise  in  the  backseat. ter  power  return.  Any  noise  in 

the  signal  will  lienee  directly  affect  the  result  for  e  ;  the  problem 

s 

will  increase  with  denser  clouds  and  longer  ranges,  due  to  the  low  SNR. 
Compounding  the  problem  is  the  uncertainty  in  q.  From  Eq .  (14),  the 
uncertainty  in  E^  is 

i-(i+*0 

(17) 
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The  behavior  of  Eq.  (17)  is  clarified  by  assuming  qz  <<  1;  then 


/  2  \  " 

Ae 

s 

i.  <1  (  “  M 

es 

q 
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s  '  s  '  J 

(18) 


Equation  (18)  shows  that  the  estimates  for  the  extinction  coefficient 
are  good  for  larger  values  of  the  smoke  transmittance.  The  second 
term  in  Eq.  (18)  is  a  major  source  of  error  only  when  Ey  is  very 
small  ('•<  Aq). 

The  quantity  E^  is  very  small  mainly  in  two  regions — in  front  of 

and  behind  the  smoke  cloud;  inside  the  cloud,  it  is  usually  much 

2 

larger  than  Aq.  In  front  of  the  smoke  cloud,  T  equals  unity,  so  the 

8 

second  term  vanishes.  Beyond  the  cloud,  theoretically  returns  to 

zero;  but  because  Efi  is  calculated  from  signals  that  have  traveled 

through  the  cloud,  the  SNR  is  very  low,  the  estimate  cf  will  suffer 

large  uncertainties,  and  the  second  term  can  become  large.  Errors  of 

a  few  orders  of  magnitude  for  e.  beyond  the  cloud  are  not  unexpected; 

fa 

but  since  e  is  near  zero,  the  absolute  effect  is  minor, 
s 
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Figure  3  plots  the  dominant  contribution  of  the  error  in  the 

estimates  of  the  extinction  coefficient  due  to  th.  uncertainty  in  q. 

Comparison  of  Figs.  2  and  3  shows  that  one  can  obtain  values  of  the 

transmittance  with  errors  less  than  the  error  in  q — i.o.,  for  values 

of  T  >  .6:  such  is  not  the  case  for  the  values  of  the  extinction  co- 
s 

efficient.  The  least  error  in  the  estimate  of  t lie  extinction  coef¬ 
ficient  is  at  best  equal  to  and  is  generally  greater  than  t lie  error 
in  q. 

> 


I 


» 


» 


» 
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TV.  OBSCURANT  CONCENTRATION  AND  CL  VALUES 


With  a  functional  value  for  the  extinction  coefficient  and  cer¬ 
tain  other  parameters  of  the  smoke  constituents  in  a  localized  region 
of  space,  it  is  theoretically  possible  to  infer  the  spatial  concentra¬ 
tion  throughout  the  region  of  nonzero  lidar  return.  To  compute  the 
concentration,  set 


and 


e  Jz) 

b 


00 

/dr  ns(r, 


2)0s(r) 


OO 

gs(z)  “  P0  f  dr  na(r*  z)  (t  rJ)  * 


(19) 


(20) 


where  n^r,  z)  is  the  number  of  smoke  particles  of  radius  r  at  range  z 
per  unit  volume,  and  0  (r)  is  the  extinction  cross  section  for  particles 
of  radius  r.  The  density  of  the  smoke  constituents  is  (in  grams  per 
cubic  centimeter),  where  the  smoke  constituents  are  assumed  approxi¬ 
mately  spherical.  The  total  number  of  smoke  particles  per  unit  volume 
N  (z)  at  range  z  can  be  written 


Ns(z) 


00 


z)  • 


Then 


e  (z)  -  N  (z)  (o  > 

b  S  b 


(21) 


(22) 


and 
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=  P0NS(z)(vs>  ’  (23) 

where  (o  )  is  the  average  smoke  extinction  cross  section  and  (v  )  the 
s  s 

smoke  particle  volume  for  the  distribution  function  n,(r,  z)  .  Lot 

k  -  P0(v8)/(as)«  s0 

C  (z)  =  ke  (z)  (24) 

gs  s 

and 


CL 

s 


C 


z 


(25) 


To  use.  the  formulas  for  C  and  CL  necessitates  estimating  k  for 

gs  s 

a  localized  region  in  the  smoke  cloud  along  the  propagation  path  of 

the  Ildar  beam,  where  values  for  e  ( z )  can  be  calculated.  There  are 

s 

two  main  methods  of  obtaining  values  for  k.  The  first  considers  the 

various  terms  for  k  separately.  For  a  known  smoke  cloud,  since  the 

material  of  the  assumed  single  scattering  component  is  known,  is 

known.  Then  from  local  measurements  of  n  (r,  zn)  and  o  (r),  (v  ) , 

(cr^),  and  k  can  be  calculated.  In  the  second  method,  the  value  of  k 

is  found  by  making  a  direct  measurement  of  C  at  some  z  and  dividing 

the  value  by  ey(z)  at  that  point. 

Accurate  values  for  C  and  CL  require  accurately  measuring 

gs  s 

n,(r,  z)  and  0  (r)  in  the  first  method.  However,  accurate  measure- 

S  8 

ments  of  n(r,  z)  seem  quite  difficult  to  obtain.  For  example,  in 
s 

the  report  from  Smoke  Week  II  [Farmer  et  al.,  1979],  comparison  of 


data  from  a  particle-sizing  interferometer  and  the  measurements  made 
in  the  Signature  Characterization  Facility  showed  "radically  differ¬ 


ent  size  distributions."  Compounding  the  errors  in  ns(r,  z)  are  those 

associated  with  uncertainties  in  p  and  e  (z). 

u  s 


A  simple  local  measurement  of  C  is  subject  to  minimal  error; 

gs 


uncertainty  in  C  g  and  CLg  at  other  ranges  will  be  due  mainly  to  the 


uncertainty  in  e^(z).  The  second  method  of  directly  measuring  the 
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mass  concentration  is  therefore  less  prone  to  uncertainty  than  the  first, 
but  it  can  be  improved  bv  measuring  e  (2)  along  with  C  (z) .  That  pro- 

25  gs 

cedure  will  obtain  a  better  value  of  C  /e  ,  as  well  as  check  the  lidar 

gs  s* 

model  by  comparing  the  calculated  with  the  locally  measured  e  .  A 

s 

further  Improvement  would  be  to  take  more  than  one  local  point  measure¬ 
ment  along  the  path  of  propagation,  then  use  the  data  to  verify  the 
model;  if  the  measurements  showed  variation  from  the  model,  the  model 
assumptions  could  be  modified. 

In  any  case,  it  is  clear  from  Eqs.  (24)  and  (25)  that  relative 
values  of  the  concentration  and  CL  are  obtained  directly  from  the 
extinction  coefficient.  The  relative  values  alone  are  useful  for 
studying  the  spatial  and  temporal  evolution  of  the  obscurant  cloud, 
even  without  an  absolute  calibration. 
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V.  APPLICATION  OF  MOOLI.  TO  DATA 


We  applied  the  technique  for  calculating  the  various  smoke,  cloud 
parameters  outlined  above  to  the  10.6  ym  bac.kscatter  data  (Fig.  1)  and 
the  transmittance  measurements  of  event  F-6  (where  artillery  shells 
were,  fired  to  produce  a  dust  cloud)  in  the  D1RT-I  program  [van  dor  Laan, 
1979  J.  The  10.6  ym  bac.kscatter  datu  at  T  “  0  show  the  ambient  air 


bnckscatter  amplitude  as  a  function  of  the  range  Pa(z),  where  the  two 


spikes  are  due  to  reflectors  for  distance  calibration  and  transmittance 

measurements.  At  T  +  19s  the  2000  m  target  return  is  not  visible.  By 

T  +  2  min,  the  target  return  is  clearly  visible,  the  cloud  is  still 

substantial,  and  there  is  a  fairly  good  signal  between  the  cloud  and 

the  target.  Since  we  had  a  particle  size  distribution  measurement  from 

event  F-6  at  approximately  V  +  2  min  [Loveland  et  al . ,  1979],  the 

T  +  2  min  backscatter  data  were  chosen  for  P , (z) ,  so  that  T,(z)  and 

e  (z)  as  well  as  C  (z)  and  CL  (z)  could  be  calculated, 
s  gt?  s 

Equation  (1)  was  used  to  obtain  a  theoretical  curve  fit  to  the 
ambient  backscatter  data  (without  the  two  spikes  from  the  reflectors) 


I’  (z).  Employing  Kq.  (1)  in  the  curve  fit,  the  explicit  expression 


for  Chz)  must  be  used,  and  we  must  as  sum-'  b  Gaussian  aperture  and 


field  of  view,  as  well  as  that  the  pulse  is  short  and  0  and  e  are 


constant.  Several  runs  were  made  to  obtain  a  curve  fit  to  P  (z).  The 

a 


best  fit  is  shown  in  Fig,  A,  with  results  for  the  parameters 


A0  =  1.8  x  10^,  c  =1.1  km  \  and  b  =  21.3  cm. 
a 


At  close  ranges,  much  of  che  backscatter  is  not  collected  by  the 
receiver  until  the  field  of  view  overlaps  the  beam.  That  is  because  of 
the  receiver/ transmitter  separation  b  and  the  associated  f ield-of-view 
and  beam  divergence.  Figure  A  shows  that  "he  bac.kscatter  power  ampli¬ 
tude  increases  until  it  peaks  (region  of  overlap  of  the  field  of  view 
with  the.  beam),  thereafter  monotonies  11  y  decreasing  with  range — as 
expected.  Figure  5  illustrates  the  beam  and  f ield-of-vluw  overlap, 


with  0  =*  1.2  mrad  and  0^  =  1.23  mrad;  the  overlap  is  at  approximately 


z  =  .175  m,  which  is  consistent  with  the  peak  shown  at  approximately 
180  in  in  Fig.  A. 


The  first  quantity  calculated  was  the  smoke  transmi ctance .  Hie 
ratio  P  (z)/P  (z)  was  known  as  a  function  of  the  range  z  from  the  back- 

C  cl 

scatter  data.  A  reasonable  value  for  a  (tv  .015)  can  be  found  from 

s 

the  results  of  Smoke  Week  II  [Sztankay  et  al.,  1979],  but  not  for  [i  . 

cl 

Method  Ub  was  therefore  applied  to  find  a  value  of  q.  A  relative 
target  return  transmission  of  R^/R^  “  .74  at  a  target  range  «  2000  m 
was  used  to  iterate  for  a  value  of  q,  using  Eq.  (13).  (R,/R  *  .74  is 

consistent  with  P  (z  )/Pa(z  )  =  .74.)  The  value  of  q  calculated  by 

C  t  U  Q 

this  method  is  7.3  x  10“'?  cm“^;  with  a  E-  .015,  3  is  5.5  x  10 

-1  -1  f!  a 
str  cm  .  Substituting  q  into  Eq .  (13),  the  smoke  transmittance 

can  be  estimated  for  all  values  of  the  range,  with  the  results  shown 

in  Fig.  6.  Figure  6  also  plots  T  for  Aq/q  <  .1  (dark  shading)  and 

s 

Aq/q  rC  .5  (diagonal  shading),  illustrating  that  the  error  in  Ty  is 
about  half  the  uncertainty  in  Aq/q.  Numerically,  AT/T  <  .05  for 
Aq/q  ?S  .1  and  AT/T  <  .3  for  Aq/q  <  .5;  the  approximation  to  the  error 
[Eq.  (16)]  gives  a  slight  underestimate. 

Once  the  smoke  transmittance  and  q  are  known,  using  Eq .  (14)  to 
calculate  c^(z),  the  local  value  of  the  smoke  extinction  coefficient 
at  runge  z,  is  quite  simple.  The  value  of  e  is  then  substituted 

9 

into  Eq,  (24)  to  derive  the  mass  concentration.  To  calculate  C  in 

gs 

absolute  units,  the  value  k  must,  first  he  found. 

The  results  of  measuring  the  particle  size  distribution  for  event 
E-G  of  DTRT-I,  conducted  at  approximately  T  =  2  min,  are  reproduced 
in  Fig.  7.  The  data  are  averages  over  a  10  sec  interval;  since  the 
velocity  of  the  probes  was  10  m/sec,  the  range  interval  is  100  m 
(assumed  to  be  in  the  center  of  the  cloud  z  ) .  Hie  particle  size 
distribution  (for  the  smoke)  n  (r,  z^)  as  estimated  from  Fig.  7,  where 
r  ranges  from  1  -  200  pm,  is 


n (r ,  z^) 


5  7  14 

3.4  X  10J  (pm) 

3.14  liter 


(26) 


Letting  pQ 
1977],  and 
in  a  value 


=  2.5  gm/cc, 
substitut ing 
for  C 


a  value  commonly  used  for  dust  [Patterson, 
that  value  and  Eq .  (26)  into  Eq.  (20)  results 
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Fig.  6--Sensitivity  of  inferred  transmittance  to  uncertainties  in  q 
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Fig.  7--Partic1e  size  distribution  measurement  for  event  F-6,  D1RT-I 
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C  (z  ) 
gs  c 


-4  3 

3.8  x  10  gm/liter  -  .38  gm/m  , 


(27) 


which  agrees  with  the  dust  results  from  Smoko  Week  II  [Farmer  et  al., 
1979].  The  above  value  of  C  is  divided  by  an  average  smoke  extinc¬ 
tion  coefficient  in  the  interval  about  z  to  obtain 

c 


3  3 

k  *  9.6  cm-gm/ liter  =  9.6  *  10  cm- gm/m 


(28) 


The  values  of  C  (z)  can  be  evaluated  for  the  full  range  0  to  2000  m. 

gs 

Figure  8  piots  the  calculations  for  the  smoke  extinction  coefficient 
(scale  on  left  side)  and  the  smoke  mass  concentration  (scale  on  right 
side),  where  scale  for  the  range  lias  been  expanded.  Also  plotted  are 
the  uncertainties  with  Aq/q  £  .1  (dark  shading)  and  Aq/q  <  .5  (diagonal 
shading).  The  *iumerical  calculations  show  that  for  Aq/q  £  .  1,  Ae  ,/e 

S  B 

is  about  .15  in  the  region  of  large  c  — very  close  to  the  theoretical 

B 

calculation  [Eq.  (18)]  of  about  .14.  The  numerical  analysis  also  con¬ 
firms  the  theoretical  prediction  for  the  behavior  of  At  /e  given 

4  S3 

small  values  of  e^.  may  he  as  great  as  10  ,  but  only  when 

<<  1;  the  plot  confirms  that  such  errors  will  not  cause  problems. 

For  values  of  Aq/q  ~  .5,  errors  in  Ac  /e  and  AC  /C  are  only  on  the 
11  s  s  gs  gs 

order  of  1. 

As  for  the  question  of  a  local  versus  a  particle  distribution 

measurement  of  C  ,  consider  the  data  for  calculating  ny(r,  z)  in 

Fig.  7.  Choosing  a  straight-line  log  fit  to  the  data  is  subjective; 

an  analysis  of  the  possible  error  shows  that  the  range  for  k  is  about 

27.8  to  1.9  cm-gm/liter— which  translates  to  possible  +190  percent  to 

-80  percent  errors  in  the  calculation  for  C  ^  before  any  errors  due 

to  Aq  are  considered,  A  direct  local  measurement  of  C  will  signifi¬ 

es 

cantly  reduce  such  errors. 

The  final  calculation  is  for  CLy  [Eq.  (25)],  where  the  values  of 

C  are  substituted  and  the  integration  is  carried  out.  The  values 
gs 

for  CL  as  a  function  of  range  are  given  in  Fig.  9,  showing  also  the 
uncertainty  due  to  Aq/q  £  .1  and  Aq/q  <  .5.  Errors  from  C  (z)  propa- 
gate  to  errors  in  CL;  added  errors  due  to  range  are  insignificant. 


Extinction  coefficient  (1/km) 


concentration  lO  ^gm/liter 
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VI.  CONCLUSIONS 


The  analytic  results  show  that  with  only  Ildar  backscatter  data, 
the  transmittance  and  the  extinction  coefficient  for  smoke  can  be 
calculated.  With  a  two-way  smoke  transmittance  of  .5,  the  estimated 
errors  are  AT  / T  rs  .5  Aq/q  and  Ae/ti  a  2Aq/q;  to  achieve  good  esti- 
mates  of  the  uncertainty  in  q  must  therefore  he  minimized.  To 
estimate  the  mass  concentration  and  CL  values  in  absolute  units  requires 
Information  on  the  obscurant  particulates.  The  usual  quantity  mea¬ 
sured  is  the  smoke  particle  size  distribution  n_(r,  z) ,  but  because 
of  the  difficulty  of  accurately  measuring  n  (r,  z)  and  the  compounded 
errors  from  other  measurements  needed  to  calculate,  k  (the  ratio  of 
mass  concentration  to  extinction  coefficient),  we  argue  in  favor  of  a 

direct  local  measurement  of  C  .A  point  measurement  of  e  would 

gs  s 

further  decrease  the  uncertainty  in  k.  Additional  local  measurements 

at  different  ranges  within  the  smoke  cloud  would  yield  information  to 

verify  the  model.  Figure  10  diagrams  the  analysis. 

Data  for  event  F-6  of  DIRT- I  give  quantitative  results  for  T  (z), 

s 

e  (z),  C  (z),  and  CL  (z).  The  results  show  maximum  errors  of  less 

a  gS  b 

than  5  percent  in  T  (z)  and  less  than  20  percent  In  c  ,(z)  for  Aq/q  *  .1. 

o  S 

For  the  same  value  of  Aq/q,  the  errors  in  C  (z)  and  CL  (z)  were  less 

g  s  a 

than  20  percent  and  12  percent,  respectively.  However,  estimates  show 

a  possible  200  percent  additional  error  due  to  the  uncertainty  in 

n  (r,  z)  and  p„,  used  to  estimate  k  in  calculating  C  and  CL  .  Such 
s  U  gs  s 

a  large  error  would  not  occur  with  a  local  measurement  of  C 

gs 

To  summarize,  the  principal  conclusions  of  the  study  are  that 

•  The  inversion  algorithm  can  produce  reliable  estimates  of 
smoke  or  dust  transmittance  and  extinction  from  all  points 
within  a  cloud  for  which  a  resolvable  signal  can  be  detected. 

•  A  single  point  calibration  measurement  can  convert  extinction 
values  to  mass  concentration  for  each  resolvable  signal  point. 

•  Having  a  target  behind  or  to  the  side  of  the  obscurant  cloud 
aids  In  normalizing  the  transmittance  and  extinction  estimates. 
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APPENDIX 


2 

From  Eq .  (2),  it  is  obvious  that  T,  is  bounded  by  0  and  1;  but 

2  ^ 

in  the  final  expression  for  [Eq.  (13)],  such  is  not  the  case.  The 
problem  arises  in  Eq .  (12),  where  all  the  terms  seem  to  take  on  arbi¬ 
trary  value  .  Physically,  P,(a),  P  (z),  3„(z),  6„,  and  a  must  be 
greater  than  or  equal  to  zero.  Those  conditions  impose  a  set  of  in¬ 
equalities  given  by 


1  + 


as 

28« 


■  .lz  J  S  !'„(*> 


Jrjju)  \ 

dz  I 


for  all  z.  If  the  Ildar  equation  [Eq.  (1.)]  accurately  describes  the 

backseatter  power  P  and  P  ,  and  if  the  assumed  linearity  between 

backseatter  and  extinction  is  valid,  the  above  inequality  will  always 

hold  for  true  data.  Therefore,  the  solution  to  Eq .  (12)  will  yield 

values  of  T2  between  zero  and  1.  For  real  data,  there  is  noise  in  P 
s 

and  P  and  uncertainty  in  the  true  values  of  0n  and  a  ,  so  a  minor 
n  "  h 

violation  of  the  inequality  is  not  sufficient  to  invalidate  the  equa¬ 
tions  or  the  assumptions. 
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